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Abstract 

The two-loop three-gluon vertex is calculated in an arbitrary covariant gauge, 
in the limit when one of the external momenta vanishes. The differential Ward- 
Slavnov-Taylor (WST) identity related to this limit is discussed, and the relevant 
results for the ghost-gluon vertex and two-point functions are obtained. Together 
with the differential WST identity, they provide another independent way for cal- 
culating the three-gluon vertex. The renormalization of the results obtained is also 
presented. 
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1 Introduction 



Jet studies are becoming increasingly precise, both as a testing ground for QCD, and as 
a background for new physics (e.g. Higgs searches). Increasing precision, among other 
things, requires knowledge of the fundamental QCD vertices to higher loops. 

The one-loop vertices have been known for quite some time. Celmaster and Gonsalves 
presented in 1979 [0 the one-loop result for the three-gluon vertex, for off-shell gluons, 
restricted to the symmetric case, p\ = p\ = V%i m an arbitrary covariant gauge. The result 
of [l[ was confirmed by Pascual and Tarrach 0. Ball and Chiu then in 1980 considered 
the general off-shell case, but restricted to the Feynman gauge |J . Later, various on-shell 
results have also been given, by Brandt and Frenkel j|] , restricted to the infrared- singular 
parts only (in an arbitrary covariant gauge), and by Nowak, Praszalowicz and Slomihski 
H, who also gave the finite parts for the case of two gluons being on-shell (in Feynman 
gauge). The most general results, valid for arbitrary values of the space-time dimension 
and the covariant-gauge parameter, have been presented in our previous paper || . Some 
results for the one-loop quark-gluon vertex (or its Abelian part which is related to the 
QED vertex) can be found in [^]. 

The present paper is devoted to a study of two-loop corrections to the three-gluon ver- 
tex in the zero- momentum limit. This limit refers to the case when one gluon has vanishing 
momentum. The remaining two momenta must then be equal and opposite, so there is 
only one dimensionful scale, p 2 . In this limit, the renormalized expressions for QCD ver- 
tices in the Feynman gauge have been presented by Braaten and Leveille ||. Information 
about Green functions is also required for calculation of certain quantities related to the 
renormalization group equations, such as the j3 function and anomalous dimensions. The 
two-loop-order contributions to these quantities were calculated in refs. || [TO], [TTJ, |P2| , 



whereas the three-loop-order results were obtained in [|13| , |14[ . Moreover, recently the 
four-loop-order expressions became available [fLq| . 

When massless quarks are considered, the scalar functions corresponding to the coef- 
ficients of different tensor structures are in the zero-momentum limit rather simple: apart 
from non-trivial coefficients, they are given by p 2 raised to some power (determined by 
the dimension of space-time). Also, the tensorial structure is considerably simpler than in 
the general case. Although the zero-momentum limit has limited physical applications, it 
serves as an important reference point, against which more general results can be checked. 

With one gluon momentum vanishing, there are two Ward- Slavnov- Taylor (WST) 
identities, one corresponding to the vanishing momentum, and one corresponding to the 
finite momentum. The identity corresponding to the vanishing momentum turns out to 
be a differential identity. In this case, the three-gluon vertex can actually be completely 
constructed from the two-point functions and the ghost-gluon vertex, with no additional 
transverse term. 

In the present paper, we realize two ways to calculate the two-loop three-gluon vertex 
in an arbitrary covariant gauge. One of them is a straightforward calculation of all 
diagrams contributing to the three-gluon vertex at this order. Another way is based on 
using the results for the ghost-gluon vertex and the two-point functions, together with 
the corresponding WST identities. The renormalized expressions are also obtained. 



2 Preliminaries 



The lowest-order gluon propagator is 

^(?«»-«^)> (2-1) 

where £ = 1 — a is the gauge parameter corresponding to a general covariant gauge, 
defined such that £ = (a = 1) is the Feynman gauge. Here and henceforth, a causal 
prescription is understood, 1/p 2 — > l/(p 2 + iO). 
The three-gluon vertex is defined as 

r ™ 3 (Pi^2,p 3 ) = -i 9 f m T m (p uP2 ,p 3 ), (2.2) 

where f aia2a3 are the totally antisymmetric colour structures corresponding to the adjoint 
representation of the gauge group (for example, SU(iV) or any other semi-simple gauge 
group). In fact, also completely symmetric colour structures d aia2a3 might be considered, 
but they do not appear in the perturbative calculation of QCD three-point vertices at the 
one- and two-loop level. Since the gluons are bosons, and since the colour structures f aia2a:i 
are antisymmetric, r AllM2Ai3 (pi,p 2 ,P3) must also be antisymmetric under any interchange 
of a pair of gluon momenta and the corresponding Lorentz indices. 

When one of the momenta is zero, the three-gluon vertex contains only two tensor 
structures^, 

P (2-3) 

In this decomposition, we basically adopt the notation of || for the scalar functions Ti(p 2 ). 
The first tensor structure on the r.h.s. of eq. ( |2.3p corresponds to the lowest-order vertex. 
There is the following correspondence between the functions Tj and the scalar functions 
A and C used in [|J (cf. also in |§): 

T^p 2 ) <- A(p 2 ,p 2 ;0), T 2 (p 2 ) <- -2p 2 C(p 2 ,p 2 -0). (2.4) 

At the lowest, "zero- loop" order, the Yang-Mills term of the QCD Lagrangian yields^ 

rf°) = l, T 2 (0) = 0. (2.5) 

For a quantity X (e.g. any of the scalar functions contributing to the propagators or 
the vertices), we shall denote the zero- loop-order contribution as (cf. eq. Q2.5D ), the 
one-loop-order contribution as X 1 - 1 ' , and the two-loop-order contribution as X^ 2 \ In this 
paper, as a rule, 

X (L) = X (L,0 + X (L, q ) ^ ( 2 g) 

where X^ L '^ denotes the contribution of gluon and ghost loops in a general covariant 
gauge (|2.1|) (in particular, X^ L '^ corresponds to the Feynman gauge, £ = 0), while X^'^ 
represents the contribution of the quark loops. 

1 This is a corollary of the differential WST identity, see in section 3. 

2 We include the contribution T^ "* = 1 into the definition of T\{p 2 ), eq. (|2.3|). 



The ghost-gluon vertex can be represented as 

rZ a2a3 (Pi,P2;Ps) = -ig r a2az Pi" r^ 3 ( Pl ,p 2 ;p 3 ), (2.7) 

where p\ is the out-ghost momentum, p 2 is the in-ghost momentum, p 3 and // 3 are the 
momentum and the Lorentz index of the gluon (all momenta are ingoing). For r m , the 
following decomposition was used in ||: 

f W3 (pi,p 2 ;P3) = 9w 3 a (P3,P2,Pl) -P3 M P2 M3 KP3,P2,Pl) +Pl M P3 M3 c(p 3 ,P2,Pl) 

+P3^Pi^d(P3,P2,Pi) +Pi M Pi M3 e(p 3 ,P2,Pi). (2.8) 

At the "zero-loop" level, 

r^] 3 = 9^31 (2-9) 

and therefore all the scalar functions involved in ( |2.8p vanish at this order, except one, 
a<°) = 1. 

We shall need the results for the ghost-gluon vertex fl2.8| ) for two different configura- 
tions: (i) when the gluon momentum, p 3 , is zero and (ii) when the in-ghost momentum, 
P2, is zero. In the former case, we get 

frtu 3 (-p,p;0) = g^a 3 (p 2 ) + p^p^e^p 2 ), a 3 (p 2 ) = a(0,p,-p), e 3 (p 2 ) = e(0,p,-p), 

(2.10) 

whereas in the latter case we obtain 

r MM3 (p,0; -p) = g^ 3 a 2 (p 2 ) + P^P^e 2 (p 2 ), a 2(p 2 ) = a(-p,0,p) } e' 2 (p 2 ) = e'(-p,0,p), 

(2.11) 

with 

e'(p 3 ,P2,Pi) = e(p 3 ,p 2 ,p 1 ) - c(p 3 ,p 2 ,Pi) - d(p 3 ,p 2 , Pl ). (2.12) 
We shall also denote 

d 2 (p 2 )=d(-p,0,p). (2.13) 

We do not need to consider r MM3 (0,p, — p) (pi = 0) because it does not enter the WST 
identities (see in section 3). Moreover, the proper ghost-gluon vertex (|2.7|) vanishes in 
this limit, for it contains p^. 

The gluon polarization operator is defined as 



K\Z(p) = [P 2 g^ -PvPv) J(p 2 ), (2.14) 



while the ghost self energy is[] 



U aia2 (p 2 ) = 6 aia2 p 2 G{p 2 ) '. (2.15) 



-1 



In the lowest-order approximation = = 1. 

3 There was a misprint in eq. (2.8) of ||: G(p 2 ) should read [G(j 
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3 WST identity in the zero- momentum limit 

In a covariant gauge, the Ward-Slavnov-Taylor (WST) identity ]16| for the three-gluon 
vertex is of the following form (see e.g. in fl7| ): 



pgr 



9^ Pi - Pl W Pl^ 



r M3^ 2 (Pl ; P3;P2; 



(P2,ps;pi)- (3.1) 



It is easy to see that the c and e functions from the ghost-gluon vertex ( |2.8| ) do not 
contribute to this identity. 

Consider what follows from ( |3.1| ) in the limit when one of the momenta vanishes. We 
should distinguish between two different cases: when the vanishing momentum is the one 
with which the three-gluon vertex is contracted, and when it is not. In the former case, 
we obtain a differential identity, whereas in the latter case we get an ordinary identity. 

In the differential case, we should consider p 3 = 5 — > 0, p\ = p, p 2 = — p — 5. We 
do not need the terms of order 5 2 and higher. In particular, G(5 2 ) = G(0) + 0(5 2 ) and, 
for massless quarks, G(0) = 1. When we expand the r.h.s. of eq. (|3.1|) in 5, the lowest 
("constant") term disappears, so only the term linear in 5 is relevant. Differentiating both 
sides with respect to <5 M3 and putting 5 = 0, we get 



H1H2H3 

(p, -p, 0) = (2g 

IJ.1IJ.2PH3 9fJ.lfJ,3Pll2 9lJ,2lJ.3PHl) K^P ) 



P 2 d 2 (p 2 ) AP 2 )G(0) 



+2p M3 \9h 



1H2 



P\JlPjJ2 
P 2 , 



p 2 d 2 (p 2 )+a 2 (p 2 )-p 



,da 2 (p 2 
dp 2 



J(P 2 



2 / 2^ J (P 2 ) 
-P W ) j 2 

dp 



G(0), 



(3.2) 



where the functions a 2 (p 2 ) and d 2 {p 2 ) are defined in eqs. ( |2.11| ) and ( |2 . 1 3| ) , respectively. 
The function a 2 (p 2 ) is defined as 



d 

a 2 {p 2 ) = p la a(p 3 ,-px -p 3 ,pi) 

■ rlfT Pl=~P3=P 

It can be calculated directly at the diagrammatic level (see in section 5). 

Considering contraction with a non-zero momentum, we get from eq. ( |3.1| ) 



(3.3) 



p 



mr MiM2M3fe -P.0) = -J{p 2 )G{p 2 )a 3 {p 2 ) (g^ 3 p 2 - PnPn 



(3.4) 



where a 3 (p 2 ) is defined in eq. (|2.10| ). Contracting eq. (|3.2|) with p Ml we get a different 
representation which should be equal to the r.h.s. of eq. ( |3.4p . Therefore, the following 
relation should hold: 



G(0) a 2 {p 2 )-p 2 d 2 {p 2 ) =G(p 2 )a 3 (p 2 ). 



(3.5) 



Using eq. (|3.5|), the differential WST identity (|3.2|) can be re- written in a way which 
involves just the a functions from the ghost-gluon vertex: 



HlH2H3 

(p, -p, 0) 



Phi I 9h2H3 



+ 2 P M3 9hi 



H2 



PjmPjjv 

P 2 



G(0) 



a 2 {p 2 



PH2PH3 



p 
d 

dp 2 



+ Ph2 9h 



1H3 



(p 2 J(p 2 ))-p 2 J(p 



PuiPm \ 

p 2 J 

2 da 2 (p 2 



a,{p 2 )G{p 2 )J{p 2 



dp 2 



+ a 2 {p 2 )J{p 2 



. (3.6) 



For the scalar functions Ti(p 2 ), the WST identity gives 



Up 2 ) = a 3 (p 2 ) G(p 2 ) J(p 2 



T 2 (p 2 )=2T 1 (p 2 )-2G(0) 



p'Ap^^^ + MpW) 



(3.7) 
(3.8) 



Therefore, the differential WST identity makes it possible to define the whole three- 
gluon vertex (not only its longitudinal part) in terms of two-point functions and the 
ghost-gluon vertex. Moreover, it can be used as another independent way, in addition to 
the direct calculation, to obtain results for the three-gluon vertex. 



4 Results for the three-gluon vertex 

We shall use dimensional regularization fl8 |, with the space-time dimension n = 4 — 2e. 
The results for unrenormalized one-loop contributions to the scalar functions T\{p 2 ) and 
T 2 (p 2 ) (in arbitrary space-time dimension) can be found in ref. ||, eqs. (4.30), (4.31), 
(4.33) and (4.34). Expanding them in e we getf] 



107 



27 



£-^ 2 )} + 0(e 2 ), (4.1) 



<i«)fj\ - rr 9 2 V ,_ J 4 ^ 20 112 i 



= (-^~ e H - 2 ^ + 1^ 2 + £ (- f - e + 1^ 2 ) } + o(e2) ' (4 - 3) 



(? 2 ?7 . 2W [8 , 40 



*TV) = Tj^r. (~P 2 r {I + f e} + 0(. 2 ). (4.4) 



In these equations, we use the standard notation Ca for the eigenvalue of the quadratic 
Casimir operator in the adjoint representation, 

f-acdfbcd = Ca S ab = N for grQup )_ (4_ 5 ) 

Furthermore, 

T = N f T R , T R = | Tr(I) = \ , (4.6) 

where I is the "unity" in the space of Dirac matrices (we assume that Tr(J) = 4), Nf is 
the number of quarks and 

" - 3?^ r(3 - *> = w^T) r(1 + E) = (' - + °^>) ■ < 4 - 7 > 



4 In all unrenormalized expressions given in sections 4-7 and in Appendix A, the bare quantities g 2 = g 2 B 
and £ = £b are understood, i.e. the same as those given in the lowest-order functions 
the renormalization is discussed, these bare quantities get a subscript "B" (see in section 



2j)-(M)- When 



Here 7 ~ 0.57721566... is the Euler constant. The e terms in the expressions QH|) — Q4.4|) 
are needed when these expressions are multiplied by terms which diverge like 1/e, e.g., 
for the calculation of reducible unrenormalized two-loop-order contributions. The e terms 
are also necessary for getting the renormalized two- loop-order results, see section 8. 

The diagrams contributing to the three-gluon vertex at the two-loop level are shown 
in Fig. l[]. Each diagram should be considered with two other "rotations" , corresponding 
to permutations of the external legs. The grey blob corresponds to a sum of all one-loop 
contributions to the gluon polarization operator, including the gluon, ghost and quark 
loops insertions^, cf. Fig. 2a of 0. Note that non-planar graphs do not contribute to the 
two-loop vertex, since their over-all colour factors vanish, due to the Jacobi identity (cf. 



Fig. 6 of ref. [20] where this is explained). 

When one external momentum vanishes, technically the problem reduces to the cal- 
culation of two-point two-loop Feynman integrals. To calculate the occurring integrals 



with higher powers of the propagators, the integration-by-parts procedure plj has been 
used. For the integrals with numerators, some other known algorithms [f2TJ (see also in 
p2| ) were employed. Straightforward calculation of the sum of all these contributions^ 
yields the following results for the unrenormalized scalar functions: 



ftO,_*. .^^V, 13 ! 15^1/ 311 n 29 tf2 7 



k 3 



W) -dffiw^ 

6965 1 509^ 15 _ 115. 13, 1 A , . , , 
"IBS " I C3 " 576* + 8^ - 144^ + l/ + 16*} + ®> ^ 



T (2,,) (p 2 



(4tt)' 



e U2 3^ 3 s 
1675 



72 



22 



ff 4 V 2 
(4?r) n 



-P 



2\-2e 



{M-i6Cs} + (,) 



(4.9) 



2\-2e 



(47T) n 



22 11. 8^ 2 7^o 

£ + -£ £ 

3 6 3 16 



1013 13 1 83 

~3g- - Cs + y4 - 2?Cs - ^> 



e 2 + 7e 3 -^ 4 +0(e), (4.10) 



T M fwVrr A 2 / fl (32 16 2 \ 289 133 4 



,4 ^2 



-8C F T^(-p 2 )- 2£ + 0(^), (4.H) 



5 To produce the figures, the AXODRAW package Q was used. 

6 Here and henceforth, we do not show contributions involving tadpole-like insertions which vanish in 
the framework of dimensional regularization Jlq] . 

7 For this calculation, two independent computer programs written in REDUCE |23| and FORM Q 
were used. 



7 



where £3 = ((3) = YljLiJ 3 — 1.2020569... is the value of Riemann's zeta function; Cf is 
the eigenvalue of the quadratic Casimir operator in the fundamental representation. For 
the SU(iV) group, C F = (N 2 - 1)/(2JV). 



5 Results for the ghost-gluon vertex 

In order to check the WST identity, we need results for the ghost-gluon vertex in two 
limits corresponding to eqs. (12. lOf) and (|2.11 ). We shall also need the derivative 52 (p 2 ), 



eq. (PD. 

The relevant one-loop results (for an arbitrary n) are listed in Appendix A. Expanding 
them in e we get 

4' V) = C, j0j- 2 (-P 2 )- (1 - + \ + e) + 0(, 2 ), (5.1) 
<4V) = C, ^5 (1 " + jf + 5^} + 0(e\ (5.2) 



(5.3) 

pM' V) = Cx ^ ("P 2 )"' {i + ^ + e} + 0( £ 2 ), (5.4) 



P V< V) = CU (V)-(l - 0(2 - { j + \c) + O(e'). (5.5) 

Two-loop contributions to the ghost-gluon vertex are shown in Fig. 2. As in the case 
of the three-gluon vertex (cf. Fig. 1), non-planar graphs do not contribute (cf. ref. ||20|| ). 
Straightforward calculation gives the following results: 



«™ - ^ < vr* {^(§-^+ + j lf - ^ + 

4 2 ' *V) = (-P 2 )-* + °(<0. I 5 - 7 ) 

^W=^(V)-{Kf4f-I^) + f + S<«-i5f + >-^K)- 

(5.8) 

p 2 e?'V) = CaTj^- (- P r 2£ {~ - 4} + 0(e), (5.9) 



(4tt 



227 
+ 72" 



?"V) = C A T$t (V)- 2 ' (1 - ?) 2 {-£ - y } + 0(e), (5.11) 



^«(P 2 ) = ^(-/)-(l-fl{iG_^ + ^ 



+ S + ^-I«-^ + ^ 2+ ^ 3 } + °( £ »' < 5 - i2 » 



p Vf ' V) = C A T0^ (-p 2 )- 2e (1 - 2 [y e + y } + 0(e). (5.13) 



The derivative (|3.3j ) has been calculated in the following way. The momenta p\ and p% 
are considered as independent variables, whereas p 2 = —pi—p3- Therefore, the momentum 
Pi flows from the in-ghost leg to the out-ghost leg. An unambiguous p\ path inside the 
diagram can be chosen as the one coinciding with the ghost line. This is convenient, 
since all we need to differentiate are just two types of objects: ghost propagators and 
ghost-gluon vertices occurring along this path. In this way, we avoid differentiating gluon 
propagators and three-gluon vertices. We also avoid getting third powers of propagators. 

Technically, this was realized as follows. The list of diagrams contributing to the 
ghost-gluon vertex, Fig. 2, was taken. Then, the propagators and vertices along the ghost 
path were "marked" by introducing an extra argument (say, z). Of course, the closed 
ghost loops should not be marked. Then, the derivative with respect to z was considered, 
and the rules for differentiating the ghost-gluon vertex and the ghost propagator (with 
subsequent contraction with pi Ml ) were supplied. It is very important that we do not 
really need expressions with different momenta; we just formally differentiate along the 
ghost line, and then perform all calculations for pi = —p 3 = p, p 2 = 0. Finally, extracting 
the coefficient of g^ 3 gives the following results for the function ( |3.3| ) : 



[T>) - Ca ^Y { P} \M2 + 16 e 32^ J + £ U8 + 96 e+ 2/ 32^ 

3085 1. 1265,. 7.. 389 ^ 2 13^ 1 VI s 

+ + -Cs + f - -£Cs + £ £ £ > + 0(e), (5.14) 

288 4 s3 576 S 8^ S 288 s 16 s 32 J v ; ' v ; 



~(2,g)/ 2\ _ f rp9 4 V 2 I „25-2 £ f 1 , l ( 17 2 1 2 

239 79 „ 7 



72 36 

6 Results for the two-point functions 



Z+ge} + 0(e). (5.15) 



Before presenting the results, let us make some general remarks. According to eq. (|2.14j ), 
the gluon polarization operator is proportional to 

J(p 2 ) = 1 + Jd)(p 2 ) + J( 2 V) + . . . (6.1) 



o 



Two-loop contributions to the gluon polarization operator are shown in Fig. 3. The gluon 
propagator is proportional to 



J(P 2 



9n 



1A'2 



P^lP^2 
p2 



;i-o 



P^lP^2 
p2 



Therefore, the transverse part of the propagator is proportional to 

Jb 2 )p = i - j (1 V) - j (2 V) + [J {1) (p 2 )f + ■■■ 

According to eq. (|2.15|) , the ghost propagator is proportional to 

G(p 2 ) = 1 + G^ 1 V) + G< 2 V) + . . . 
The ghost self energy (which is inverse to the propagator) is proportional to 



(6.2) 



(6.3) 



(6.4) 



G(p 2 



1 - G^V) - G^^V) + . . . 

i-g^(p 2 )-g^\p 2 ) + \g^(p 2 : 



+ 



(6.5) 



Note that the one-loop contribution to the ghost self energy gives -G^(p 2 ). Two-loop 
contributions to the ghost self energy are shown in Fig. 4. They give — G < - 2 - )( - irred ^(p 2 ). 
According to eq. (|6.5| ), the two-loop contribution to the ghost propagator consists of two 
parts, the irreducible one and the reducible one, 



G( 2 V) = G ^ cd \p 2 ) +G^ red \p 2 



(6.6) 



where G^ {rcd \p 2 ) = [G {1 \p 2 ] 

One-loop results in arbitrary space-time dimension are available e.g. in |(J (see also 
in Appendix A). When we expand them in e and keep the terms up to the order e, we get 



j(U)(p2 



g 2 v 

( 47r )n/2 



i-pril 



-A 



31 
~9~ 



J^)(p 2 )=T 



C 



g 2 v 

( 4vr )n/2 

g 2 v 



-v 



2\-e 



(6.7) 
(6.8) 
(6.9) 



1 fl 1 

/ — n~\~- > 

( 4vr )n/2 

Calculating the sum of one-particle irreducible two-loop diagrams contributing to the 
gluon polarization operator (shown in Fig. 3), we have obtained the following unrenor- 
malized results: 

A ^2 



J^\p 2 



c 



2 9 V 



-p 2 r 2e 

14311 
432 



1 



25 



1 



! - \-2 ' 2/ ' / 



583 113, 19 , 3 , 

I £ £ +-£ 3 

72 144 s 24 s 8 



+ C 3 + i|« + 2«3-^ + ^ + ^} +O(£ ),(6.10) 



1 n 



(4vr 



1961 142 22 2 



+CfT w^ (Vr2£ {! + f - 16Cs l + ° {£) - (6 - n) 

Calculating the sum of the contributions (Fig. 4) to the ghost self energy (with a minus 
sign, cf. eq. (|6.5| )), we obtain 



503 3 t 73. 3. 9 3 



(4tt 



3-— £ + -£ 2 ^CsKOfe), (6.12) 

32 4 s 64^ 8^ 16 s r 



Note that there is no reducible part in G^ 2 '^. The reducible part of is given by the 
square of eq. ( |6.9|) , 



G («>mV) = {I (i + i{ + + K 1 + 5 ? ) +3 + ? } + ° (£) 



(47T 

Therefore, using eq. (|6.6| ) we get 



(6.14) 



(4vr 

599 3 9 3 ,, 3 

4 



4 s 64^ 8 s 16^ J K ' K ' 



32 

7 WST identity at the two-loop level 

Due to the differential WST identity, we get the representations (|3.7| ) and (|3.8|) for the 
functions Ti(p 2 ). In the massless case, all one-loop expressions are proportional to {p 2 )~ £ , 
whereas two- loop expressions contain {p 2 )~ 2s ■ Thus, the differentiations in (|3.8j ) become 
trivial. Expanding in g 2 , we get0 

t/V)=4V) + G ( V) + ^ ( V), (7.1) 



T [ 2 \p 2 ) = 4V) [g {1) (p 2 ) + J {1) (p 2 )} + GMCpV'V) 

+af\p 2 ) + G ( - 2 \p 2 ) + J^(p 2 ), (7.2) 
T^\p 2 ) = 2TfV) - 2 [(1 - e) jW(p 2 ) + (1 + e)4V) + a?V)l , (7.3) 



3 We take into account that (in the massless case) G(0) = 1. 



I I 




(p 2 ) = 2TfV) - 2 [j<V)<4V) + j<V)2?V) 

+(1 - 2e) J (2 V) + (1 + 2e)4 2) (p 2 ) + 4 2) (p 2 ) 



(7.4) 



Substituting the expressions for ghost-gluon vertex and two-point functions, we arrive 
at the same results as given in ( |4.8|) -( [4~TT|) . 

8 Renormalization 

To begin this section, we would like to explain why the zero-momentum limit of the three- 
gluon vertex, as well as the relevant limits of the ghost-gluon vertex, are infrared finite, 
i.e. we do not get any 1/e poles of infrared (on-shell) origin. The main argument is just 
power counting. 

Consider a triple vertex Vq (part of a two-loop diagram) to which are attached the 
zero-momentum external line, together with two adjacent propagators carrying the same 
loop momentum q. In the case of a scalar (say, 3 ) theory, one would get l/(q 2 ) 2 in the 
integrand, leading to an infrared divergency. However, in QCD the vertex Vq can be either 
(i) a three-gluon vertex, (ii) a ghost-gluon vertex, or (iii) a quark-gluon vertex. Effectively, 
the power of the gluon or ghost propagator in QCD is l/(g 2 ), whereas for the massless 
quark propagator we get 1/q. Therefore, the case (iii) is infrared finite, since we get only 
1/g 2 from the two quark propagators (no g-dependent factor from the vertex). In the 
cases (i) and (ii), we get l/(g 2 ) 2 from the two gluon (or ghost) propagators. However, we 
also get a momentum- dependent factor from the three-gluon (or ghost-gluon) vertex Vq, 
which cannot contain any momentum other than q (since the external momentum is zero). 
This gives in the numerator a factor which is linear in q, so that effectively the infrared 
behaviour is just 1/q 3 , i.e. we have no infrared divergency. When the zero-momentum 
line is attached to the four-gluon vertex like e.g. in diagrams (h) and (h') in Fig. 1, we 
may also get two propagators carrying the same momentum q. However, a similar power 
counting shows that there are no infrared singularities. For example, in diagrams (h) and 
(hf) an extra momentum q appears in the numerator from the one-loop self-energy-type 
insertion. This explains why all singularities in this limit are of ultraviolet origin, and 
therefore should be removed by renormalization. 

In this paper we adopt the modification of the renormalization prescription by 't Hooft 
[p7[ , corresponding to the so-called MS scheme ||28|| . In this section (and in Appendix B), 
the notations £, a, g 2 , etc. (without subscript) correspond to the renormalized (in the MS 
scheme) quantities. In previous sections (and in Appendix A), they should be understood 
as the bare quantities £g, as, g%i e ^c 

The renormalization constants relating the dimensionally-regularized one-particle- 
irreducible Green functions to the renormalized ones, 



(ren) I J Pi 




) 



lim Z v (-,a,g 2 J T ([p 2 },a Bl g 2 Bl e) 



(8.1) 



look in this scheme like 




(8.2) 
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where a = 1 — £. In eq. (|8.1| ) /i is the renormalization parameter with the dimension of 
mass. It is assumed that on the r.h.s. of eq. ( |3.1| ) the squared bare charge g 2 B and the 
bare gauge parameter a B must be substituted in terms of renormalized ones, multiplied 
by appropriate Z factors (cf. eqs. ( |8.8| ) and (PT5|) ) . 

We use the following definitions for renormalization factors: 



r i r i 3 (pi>P2.P3) = Z\ r WAt2M3 (pi,p 2 ,p 3 ), 
nK aia2 (p) = ^n^ 2 ( P ), 



(8.3) 
(8.4) 
(8.5) 

(8.6) 

where n^ 2 (p) and Il aia2 (p 2 ) are the gluon polarization operator and the ghost self en- 
ergy, respectively. For the scalar amplitudes, eqs. (|8T5|) — (|876|) mean that J{p 2 ) and G{p 2 ) 
should be renormalized by means of Z3 and Z^ 1 , respectively. Furthermore, according to 
eqs. (|8~3| ) -(|8~4 ] ) the three-gluon amplitudes (Ti and T 2 ) should be renormalized using Zi, 
whereas for the ghost-gluon functions (03, e 3 , a 2 and e' 2 ) one should use Z\. 
The WST identity requires that 



•7) 



^3 _ ^3 

If this condition is satisfied, the WST identity is valid for the renormalized quantities, 
too. 



Using (p.7|), the bare coupling constant g B can be chosen (in the MS scheme) asf] 

Ve 7N - 



47T 



2 r / 1 ^7— 2 

3 ^3 



47T 



( 7 2 z 1 2 z 3 - 3 . 



The gauge parameter a = 1 — £ is renormalized as 
ag = ^a, so that £b 
Below we shall use the following notation: 



1-^(1-0- 



h 



(4tt) 



47r : 



where 



^1 
47r' 



i.9) 



(8.10) 



The two-loop-order results for the renormalization factors have been obtained in |10], 
PI , |12] (see also in ref. [p6|). For completeness, we list the corresponding expressions in 
Appendix B. 

Using eqs. (|4.1fl -(|I~5D, ( |4.8|) - (|4TTID , (|8.3|) and ( |B.1| ), we obtain the renormalized scalar 
amplitudes appearing in the three-gluon vertex (cf. eq. fl2.3|) ), 

35 1.1 



n 



(rcn) 



1 + h 



+h 



C 2 



18 
4021 



2«" 

4" 



4^/ 9 
2317,. 15_ 



113 



1 



1 



288 4 ,J 576 

,875 , 2(1 10, 
+C A T — + 8C3 + — f - -£ 2 



-r — r + — r 

144 s 16 s 16 s 



72 



9 



9 



160 



0(/i 3 



(8.11} 



9 The factor (e 7 /(47r)) e = exp [5(7 — ln(47r))] in eq. ( j3ij| ) represents the difference between the MS and 
MS schemes (cf. also eq. (fL7|)). 
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-i(ren) 



t7 A (-i-2e + ^) + |r 



+ h 2 



641 



C A T 
287 



157 37 . 2 
1 



9 



19 



_ r j t £/■ _ ___C Z _i _C 

36 S3 18^ 2^ S3 144^ 16^ 



18 s 9 s / 

e 4 



+ 8C F T 

+ 0{h 3 ). (8.12) 



Here and henceforth, we put p 2 = — /i 2 in the renormalized expressions. In Feynman 
gauge (£ = 0), our expressions agree with eq. (B4) from ||. However, the one-loop part 
of the result for T2 in an arbitrary (non-Feynman) gauge disagrees with eq. (A10) from 

En 

The renormalized expressions for two-point functions are 



J(ren) = 1 + h 
+ h 2 



31 
"9" 



20 



3245 287 61 2 3 3 1 4 

" 144" + C3 ~"96" e + 2eC3+ 72 e "I6 e + 16 C 



451 



10 



10 



+C A T — + 8C3 + — i - — e + C F T — - 16C: 



36 



9 



55 



+ 0(/i 3 ), (8.13) 



G (ren) = l + h( J A + tf 



/ 997 _ 3 _ 41 3 2 _ _3_ 2 
V 96 4^ 3 64^ + 8^ 16^ 3 



95 
24' 



+C(/i 3 ). (8.14) 



In Feynman gauge, eq. (|8.13|) gives the same as the first of eqs. (B3) in ref. ||. Taking 
into account that 

G" 1 ] (r ° n) = 2 - G (ron) + h 2 C\ + 0(h 3 ), (8.15) 

we have also confirmed the second of eqs. (B3) in ||, i.e. the result for the ghost self 
energy in Feynman gauge. 

The renormalized expressions for the scalar functions occurring in the ghost-gluon 

vertex are 



(ren) 



i + -hc A (i-o 



c 



2 /137 1 



V 48 



299 
"96" 



1 



16 



1 



+ 0(/i 3 ),(8.16) 



P 2 et l) = - A hC A (2 + + h 2 



^ ( T + 8 C3_ 12 e+ 16^ 3 "l6^ 



-C A T 



(8.17) 



C 



/167 
A v72" 



43 
144 



f £ 

16 s 



16 



9 



CUT(l-£) 



+ G(/i 3 ),(8.18) 



2 /(ren) 

p e\ 



-/iCa (1-0(2-0 



29 5 ^ 5 „ 

h-C3 £" 

36 8 s 36 s 



16 



16 



16 



9 



+ C(/i 3 ). (8.19) 



10 Cf. footnote 19 on p. 4101 of ||. In our notation, in the HCa part of ( |8.12| ) the term i£ 2 is missing 
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We note that these functions are in the following correspondence with the functions 
G lt2 (p 2 ) used in §, eq. (A3): 

a 3 +P 2 e-s «-> 1 + G 2 , a 2 +p 2 e' 2 ^l + G 1 . (8.20) 

Using this connection, we have confirmed the two-loop-order results for G% and G 2 in the 
Feynman gauge, eq. (B5) of ref. ||, as well as the one-loop-order results for G\ and G 2 
in an arbitrary covariant gauge, eq. (All) of ||. 



9 Conclusion 

In the limit when one of the gluon momenta vanishes, we have calculated the two-loop 
contributions to the three-gluon vertex, in an arbitrary covariant gauge. In fact, we 
needed to calculate two scalar functions, Ti{p 2 ) and T 2 (p 2 ), associated with different tensor 



structures, cf. eq. (|2.3|). Two independent ways of calculating these scalar functions have 
been realized. One of them is based on the straightforward calculation of all diagrams 
contributing to the two-loop three-gluon vertex shown in Fig. 1. 

Another way of determining T\{p 2 ) and T 2 (p 2 ) is based on exploiting the differential 



WST identity (3^). In this way, we obtain representations of the scalar functions T\{p ) 
and T 2 (p 2 ), eqs. (|3.7| ) and ( |3.8|) , in terms of the functions occurring in the ghost-gluon 



vertex (Fig. 2), its derivative (|3.3j), the gluon polarization operator (Fig. 3) and the ghost 
propagator (cf. Fig. 4). We have calculated all these functions and confirmed the result 
of the straightforward calculation. 

The construction of the differential WST identity is of a certain interest, since in 
this limit it completely defines the three-gluon vertex, without leaving any "undetected" 
transverse contributions. 

We have constructed renormalized expressions for all Green functions involved. Note 
that in the zero-momentum limit the three-gluon vertex has no infrared (on-shell) singu- 
larities, this is a "pure" case for performing the ultraviolet renormalization. 

The obtained results can be considered as the first step in constructing expressions 
for the QCD vertices in more complicated cases, including on-shell configurations and the 
general off-shell case. In principle, the techniques for calculating the corresponding scalar 
integrals are already available p{| ^0] . 
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Appendix A: One-loop expressions for arbitrary n 

At the zero-loop level, we have 

a (°) = 4°) = 1) s (o) =Q) rf (o) = Q) j{ o) = G (o ) = 1) (A _ i; 
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and the r.h.s. of eq. (|3.2|) restores the zero-loop result for the three-gluon vertex, 

^ 'mi/*ja*3 ~Pi 0) = ^5 'm^P V3 ~~ 9tnn3Pfi2 ~ 9n2tJ.sPtJ.i- (A-2) 

At the one-loop level, the expressions obtained in || give the following results in the 
zero-momentum limit: 

9 2 V C A 



4V) = 

p 2 4V) = 

4V) = 
p 2 4V) = 

pV?V) = 

9 2 V C A 
" (4vr) ri / 2 32 

In these equations. 

<P 2 ) = 



(4tt)"/ 2 4 

( 4vr )n/2 



k( p 2 ) (n-2)(l-£), 
A n(p 2 )(n-A) [2 + (n-3)£], 



( 4vr )n/2 



^ (1 - [4(n - 3) - (n - 4)£] , 

^ k(jo 2 ) [2(n - 6) - (5n - 18)^ + {n - 4)£ 2 
^*«(p 2 )(l-0(2-0(n-4), 



(4 7r )n/2 

g 2 v ('\ 

(4 vr )n/2 

CU , 2 



(A.3) 
(A.4) 
(A.5) 
(A.6) 
(A.7) 



4V) 



(p 2 ) {8(n 2 - 6n + 10) - 2£(3n 2 - 26n + 52) + £ 2 (n - 4)(n - 6)} . 

(A.8) 



(V) 



2\(n-4)/2 



1 



(n-3)(n-4) v ^ ' e(l-2e) 
The results for two-point functions are (cf. e.g. in pS], 



:v) 



2\-e 



(A.9) 



J {l \p 2 ) 



g 2 7] k(p 2 ) 



(4tt)"/ 2 (n-1) 



4(3n-2) + 4(n-l)(2n-7)f - (n- l)(7^-4)£ 2 " 
+2T(n-2)), (A.10) 



P 



^ k(^) [2 + (n - 3)£] . 



(4^)n/2 4 

Taking into account that 

■(a 3 -p a d a )jl (1) = 4 1) -p a 4 1) +J< 1 J, 



(A.ll) 



(A.12) 



p 2 d 2 + a 2 -p 



2 da 2 
dp 2 



J + p 2 a 2 



dJ 
dp 2 



(i) 



.2j(l) ,^(1) ^2 dfl 2 



(1) 



,dJ« 



' 1 "i^i Zj ^ 7"*) 



dp 2 



= P »«*'+«f>-^«f> + ^jW(A.13) 

we have checked that eq. ( |3.2[ ) is satisfied at the one-loop level, for an arbitrary n. Fur- 
thermore, 



4V)-P 2 4V) = 4V) + ^V) 



9 2 V C A l , 



«(p ) (n-0- 



(4 7r )n/2 4 

Therefore, eq. (|3.5|) (which follows from eq. (3.4)) is satisfied at the one- loop level. 



(A. 14) 
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Appendix B: Renormalization factors 



The expressions for the relevant two-loop-order renormalization factors have been pre- 
sented in refs. |L0|, [11], (cf- also in p6| ). For completeness, we present the corresponding 
expressions hereQ: 



#1=1 + 



h 



-*T 
3 

13 7 
~8 ~ 16 



+ /i 2 C A T 



e+Sn-r 



32 



1 /71 45 



25 
12e 



^ = 1-^(1-0 + ^1(1-0 



£ V48 32" 16 
3 1 



W5_M + 1 



8 16 



-C F T 

-0(h 3 ), (B.l) 

+ 0(/i 3 ), (B.2) 



£ 



h 2 \C A T 
5 _ 1 



-It 

3 

rl + 24^ + 4^ + ilY + T6^^ 



V3 3 V 2e 
1 /23 15 



-C fT 

e 



+ 0(h 3 ), (B.3) 



1 + K 2 



h 2 \C\ 



16 



2e 2 



49 1 
48 ~ 32^ 



-1 

12e/ 



+ 0(h 3 ), (B.4) 



where e = (4— n)/2 and ft, = ^ 2 /(47r) 2 . One can check that eqs. (|B.1|) - ([B.4| ) obey the WST 
identity (|8.7| ), so only three of them are independent. Using the results for unrenormalized 
Green functions, we have performed an independent check on these Z factors^. 

The results for these renormalization factors (without fermionic contributions, i.e. 
for the pure Yang-Mills theory) were first presented in |10| (Feynman gauge) and |ll] 
(an arbitrary covariant gauge). The complete results in an arbitrary covariant gauge, 
including the fermionic contributions, were presented in [12] (cf 



also in §2§). In [0. 

the renormalization factors Z 3 and Z 3 were denoted as Z 2 and Z 2 . There was an obvious 
misprint in the last term of the expression for Z 2 where \T 2 should read y-tiV (in their 



notation, T 



Cp, Co 



Ca, tN 



■i?, U2 ^ oa, tiv T). We note that this misprint was copied over 
to the review [31 and the textbook [p5|l . In |25| , in the end of the first line of eq. (C.6) 
for Z 3 , the term a R Cp should read CqTrNj {olr is the renormalized gauge parameter, 
Cq <-> Ca). Then, in the beginning of the last line of eq. (C.5) for Z 3 , \Cq should 

read ^Cq. 



There are several misprints in eq. (2.30b) of [P. The term 2ffi (i) 



/v- 



2N 



should read y (|) | (qq is the renormalized gauge parameter, 



n <-> N f , f <-> T). In 

the previous term, ~ should read In the term involving ^, the "factor" | with the 
following bracket should be removed. In the one- loop-order part, ^ should read 



n As in section 8, the renormalized quantities £ = 1 — a, g 2 , etc. are understood. 

12 Note that the two-loop results for Z factors in the MS scheme are of the same form as in the MS 
scheme; the only difference is that g 2 in the definition of h should be understood as the renormalized 
squared charge in the MS scheme. 
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cf. eq. (2.30a). Finally, in eq. (2.31b) for Z\, the one-loop-order contribution should be 
multiplied by \, cf. eq. (2.31a). 

Using the 1/e term of the renormalization factor Z? (cf. eq. ( |8.2|) ), one can obtain the 
corresponding anomalous dimension 7r via 







7r 



(B.5) 



We have checked that in the Feynman gauge £ = (a = 1) the results for the anomalous 



dimensions 71, 73 and 73 coincide (in the two- loop approximation) with those from [|l3 



The anomalous dimension 71 is related to the others via 71 — 73 = 71 — 73 (this follows 
from the WST identity (|8.7|) and the definition ( p.5| )). Moreover, since (cf. in [p~3| ) 



27i (a, 9 2 ) ~ 7s (a, 9 2 ) ~ 27s (a, 



(B.6) 



we obtain the same result for the two-loop (3 function as those given in || |Tl|, |l2fP|. 
namely 



h 



~Y Ca + 3 T 



+ h 2 



~ c2 a + y CaT + ACfT 



+ O (V) . (B.7) 



Higher terms of the (3 function are available in refs. [13, 14, 15 
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Figure 1: Two-loop three-gluon vertex diagrams. 
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Figure 2: Two-loop ghost-gluon vertex diagrams. 
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Figure 4: Two-loop ghost self-energy diagrams. 
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